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ABSTRACT 

We have extracted over 400 clusters, covering more than 2 decades in mass, from 
three simulations of the rCDM cosmology. This represents the largest, uniform cat- 
alogue of simulated clusters ever produced. The clusters exhibit a wide variety of 
density-profiles. Only a minority are well- fit in their outer regions by the widely used 
density profile of Navarro, Frenk & White (1977, hereafter NFW). Others have steeper 
outer density profiles, show sharp breaks in their density profiles, or have significant 
substructure. If we force a fit to the NFW profile, then the best-fit concentrations 
decline with increasing mass, but this is driven primarily by an increase in substruc- 
ture as one moves to higher masses. The measured temperature-mass relations for 
properties measured within a sphere enclosing a fixed overdensity all follow the self- 
similar form, T oc M^/"^, however the normalisation is lower than in observed clusters. 
The temperature-mass relations for properties measured within a fixed physical radius 
are significantly steeper then this. Both can be accurately predicted using the NFW 
model. 

Key words: galaxies: cluster: general - cosmology. 



1 INTRODUCTION 

Clusters of galaxies are used to constrain cosmological pa- 
rameters as they are the largest gravitationally bound sys- 
tems and preserve imprints of the evolution of the universe. 
One of the most reliable predictors for the mass-function of 
collapsed objects is that of Press & Schechter (1974). Al- 
though based on naive assumptions, it has been widely used 
because of its simple form and, more importantly, because 
the mass function is in excellent agreement with many N- 
body simulations (e.g. Lacey & Cole 1993; but see Gross 
et al. 1998; Governato et al. 1999; Jenkins et al. 2000). 

In practice, cluster masses are not easily measured. For 
this reason, other properties such as X-ray luminosity or 
temperature are used to estimate the mass. The former of 
these is notoriously unreliable as it depends strongly upon 
emission from the core of the cluster; hence we will concen- 
trate on the latter in this paper. A simple application of the 
virial theorem suggests that the mean temperature of an ob- 
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ject that has collapsed to the virial radius, rvir, follows the 
scaling law 



(1) 



where Mvir is the mass within rvir and 2form is the formation 
redshift. In practice, clusters are assumed to have formed at 
the redshift that they are observed, and the virial radius is 
defined as the radius of a sphere enclosing an overdensity 
of 180 (in critical-density cosmologies) relative to the mean 
density at that redshift. 

Equation ^ is often assumed to be a rigorous theoreti- 
cal prediction, but this is not the case. It assumes that the 
constant of proportionality in Equation ^ is the same for 
low- and high-mass clusters, but there is no reason why this 
need be true. For a power-law density fluctuation spectrum 
of dark-matter particles, the theory of self-similarity tells us 
that the population of clusters at one redshift is a scaled 
version of the population at another redshift, but that is 
not the same as saying that coeval low- and high-mass clus- 
ters are similar in form. Even if the clusters do have similar 
morphologies, we may still get different proportionality con- 
stants if we measure their properties within the virial radius. 
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as this will correspond to different multiples of the half-mass 
radius in each case. 

Another deficiency of Equation ^ is that it predicts only 
the virial temperature — there is empirical evidence that the 
observed (emission-weighted) X-ray temperature of the in- 
tracluster medium, Tx, is greater than the virial tempera- 
ture (e.g. Edge & Stewart 1991a,b; BahcaU & Lubin 1994). 
Were these to be in constant ratio, then the scaling relations 
would be preserved, however this is unlikely to be the case. 
Physical processes such as shock-heating and radiative cool- 
ing act only on the gaseous component of clusters and not 
on the dark matter. Furthermore, these two processes have 
different dependencies upon density and so will break the 
expected self-similarity. 

Hjorth, Oukbir & van Kampen (1998), in a sample of 
8 clusters whose masses are inferred from weak and strong 
gravitational lensing and temperatures are measured from 
ASCA, confirm that (very roughly) Tx oc M/r. Most ob- 
servers, however, report the direct relation between Tx and 
M without an explicit dependence upon r and so that is 
what we will do here. The observations are discussed fur- 
ther in later sections when we can compare them with our 
own results. 

The scaling relation of Equation ^ has been tested using 
hydrodynamical simulations with varying deg rees of success. 



We will discuss these further in Section 6.1. On the basis 
of this, various authors have sought to constrain cosmologi- 
cal parameters from the observed X-ray cluster temperature 
and luminosity functions, for example, Lilje (1992), Viana 
& Liddle (1996), Eke, Cole & Frenk (1996), Henry (1997), 
Markevitch (1998), Eke et al. (1998), and Viana & Liddle 
(1999). 

Because the scaling laws of galaxy clusters are used so 
extensively as the basis of cosmological models, we are un- 
dertaking a series of N-body, hydrodynamical simulations 
to investigate the X-ray properties of clusters of galaxies 
in detail. It is intended eventually to produce catalogues of 
clusters covering at least two decades in mass, for several 
different variants of the CDM cosmology, and with a variety 
of physical processes acting on the intracluster medium. In 
this current paper we report on the first of these catalogues 
for a non-radiative simulation in the rCDM cosmology. The 
catalogue is generated by combining three simulations (128'^ 
particles each of gas and dark-matter) in which we resolve 
over 400 clusters at moderate resolution. Thus we have the 
correct boundary conditions for large-scale structure and a 
reasonable dynamic range with which to test the scaling re- 
lations. 

We look at the circular velocity profiles of the clus- 
ters to see whether they all have the same form, such as 
the commonly-used NEW model. We find both that indi- 
vidual clusters show a wide range of profiles and that the 
parameters describing an average cluster are a function of 
mass. Thus, there is no a priori reason to suppose that the 
cluster population should follow the scaling law described 
in Equation ^. Nevertheless, it turns out, somewhat fortu- 
itously, that the deviations of the mean cluster population 
from the scaling law are small, although the scatter from in- 
dividual clusters can be quite large. When measured within 
a fixed physical radius, however, the dependence of temper- 
ature upon mass is steeper than that given by Equation |l]. 

Because this current simulation does not include ra- 



Table 1. Run parameters for each of the simulations: 
box size/?i~^Mpc; softening//i~^kpc; dark-matter particle 
mass/h~-'MQ; minimum resolved cluster mass/?t~-'^MQ; mini- 
mum ratio of the 2-body relaxation time in the core of the clusters 
to the age of the Universe. The runs used in the catalogues are 
shown in normal font; the comparison test runs are shown in ital- 
ics. 
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68 


4.6 X 10" 


2.37 X lO^-* 


3.4 



diative cooling, it cannot be an accurate reproduction of 
the real Universe. Nevertheless, we have deliberately kept 
the model simple in order to investigate how accurately, or 
otherwise, this simple model follows the self-similar scal- 
ing laws between temperature and mas^. Future papers 
will then look at the additional effects of cooling and non- 
gravitational heating. 

We briefiy describe the simulations and cluster identifi- 



cation in Section t 
length. In Section 



and also discuss our choice of softening 
, we look at the density profiles of the 
clusters and compare them with the NEW model. We de- 
rive predictions for the scaling relations in the NEW model 
in Section ^, then compare these with those of our simulated 
clusters in Section |^. We compare our results with previous 
work in Section ^ and summarise our conclusions in Sec- 
tion 1^ 



2 NUMERICAL METHOD 
2.1 The simulations 

We have carried out three simulations with 128^ particles 
each of gas and dark matter. The cosmological parameters 
were as follows: density parameter, Q = 1; cosmological con- 
stant, A = 0; power spectrum shape parameter, F — 0.21; 
and a linearly-extrapolated root-mean-square dispersion of 
the density fiuctuations on a scale 8 Mpc, erg = 0.60. 
The Hubble parameter is irrelevant as we did not allow the 
gas to cool; the baryon fraction was set to a low value, 
Qb = 0.06 so that the gas makes only a minor contribution 
to the gravitational potential. The three simulations had dif- 
ferent box-sizes, corresponding to different mass-resolutions, 
as listed in Table |^. We also carried out two further simula- 
tions of the middle-sized box to test the effect of changing 
the gravitational softening: these are shown in italics in the 
Table. 

We use a parallel version of the Hydra N-body/SPH 
code (Couchman, Thomas & Pearce, 1995; Pearce & Couch- 
man, 1997). The simulations were executed on the Cray T3E 
at the Edinburgh Parallel Computing Centre as part of the 
Virgo Consortium's programme of investigations into large- 
scale structure. 



t Because we are complete in mass, we quote T-M rather than 
M-T relations 
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2.2 Cluster identification 

Initially we identify clusters in our simulation by searching 
for groups of dark matter particles within an isodensity con- 
tour of 180, as described in Thomas et al. (1998). We work 
with a preliminary catalogue of all objects with more than 
30 particles, then retain only those which have a total mass 
within the virial radius exceeding Mum, corresponding to 
500 particles of each species. The use of a small mass for 
the preliminary cluster selection ensures that our catalogue 
is complete. We have checked that using a different isoden- 
sity threshold, a different selection algorithm, or using gas 
particles instead of dark-matter particles to define the clus- 
ter, leads to an almost identical cluster catalogue — the only 
difference being the merger or otherwise of a small number 
of binary clusters. 

We define the centre of the cluster to be the position of 
the densest gas particle. This will usually correspond to the 
peak of the X-ray emission and has the advantage that it is 
independent of the cluster selection method. 



2.3 Substructure statistic 

Much of the modelling that we will do on the clusters sup- 
poses that they are smooth and spherically-symmetric. In 
practice most clusters show some degree of substructure. 
Following Crone et al. (1996) and Thomas et al. (1998), 
we measure substructure by comparing the positions of the 
density maximum, rd, and the centroid, of the cluster, 
where the latter is averaged over all particles within an iso- 
density contour of 180 times the background density. More 
specifically 



5" = 



^180 



(2) 



where rigo is the radius of a sphere enclosing a mean density 
equal to 180 times the background density. 

In some of the figures that follow, we indicate the clus- 
ters with the more prominent substructure, S > 0.2, by 
plotting them with open symbols: 14 per cent of the clusters 
fall into this category. 

2.4 Choice of softening 

In any N-body simulation, it is necessary to introduce a 
gravitational softening in order that 2-body interactions do 
not become important. Thomas & Couchman (1992) esti- 
mated the minimum ratio of the 2-body relaxation time to 
the age of the Universe (which occurs in the core of the 
cluster, close to the softening radius) to be 



to 



1.3 N,' 



soft 10ft~^Mpc 
10/i~ikpc box 



(3) 



where we have put in parameters appropriate to the current 
simulations. Here Ns is the number of particles within ^/e 
times the softening length. The values given in Table |l] are 
for the smallest clusters extracted from each run; the relax- 
ation time scales as T^^^ for larger clusters because these 
have more particles within the softening length. 

In a similar calculation to that of Thomas & Couchman, 
Steinmetz & White (1997) estimated the mean 2-body heat- 
ing timescale for the cluster as a whole. They found that 



q Q 



+ +++ + ++ + 



Figure 1. The ratio of the mean specific energy (thermal plus ki- 
netic) of all particles to the gas specific energy within the virial ra- 
dius for three different values of the softening: crosses, 10/i^^kpc; 
circles, 50/i~^kpc; stars, 100/i~^kpc. 



(tr) 



3 3N 



TV 32 In A 



(4) 



where N/2 is the number of particles within the half-mass 
radius, Rh, of the halo, InA « 3 is the Coulomb logarithm, 
and tc = 2TiRh/vc is the orbital period at Rh, where Vc is the 
circular velocity. The dependence upon the softening is much 
smaller than in Equation ^ which was measuring numerical 
relaxation in the core of the cluster. Putting in numbers 
appropriate to our halos, we find the heating timescales for 
the lowest-mass clusters in each run are approximately 10 
times the age of the Universe (the heating rate scales roughly 
in inverse proportion to mass for higher-mass clusters) . Thus 
we might expect the gas to be heated by up to 10 percent 
relative to the dark matter (the results of the tests, reported 
below, suggest that the heating is slightly smaller than this). 

Because theoretical estimates of the numerical heat- 
ing are so uncertain, we tested the sensitivity of our re- 
sults to changes in the softening. Three simulations of the 
112.9 /i"^Mpc box were carried out which differed only in 
having softenings of 10, 50 and 100 /i~^kpc, as shown in Ta- 
ble ^. For each cluster detected in the simulations, we mea- 
sured the mean specific energy of particles within the virial 
radius (the radius of a sphere that encloses 180 times the 
mean density). Figure |l| shows /3i8o, the ratio of the mean 
specific energy of all particles to that of the gas particles. 

As can be seen from the Figure, the effect of using 
too small a softening can be severe. Putting soft=10 /i~^kpc 
into Equation ^ gives tr,min/io ~0. 12-0.24 for the lowest-to- 
highest mass clusters shown in the Figure. The Steinmetz 
& White criterion makes only a very slight distinction be- 
tween the three runs with different softening (which enters 
Equation ^ only via the Coulomb logarithm) and so should 
be regarded as a necessary but not sufficient condition to 
prevent artificial heating. 

Figure ^ shows the profile of /3 as a function of radius 
(i.e. /3(r) is the ratio of the mean specific energy of all parti- 
cles to the mean specific energy of gas particles within spher- 
ical shells about the cluster centre) . Because individual pro- 
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Figure 2. The ratio of the mean specific energy (thermal plus 
kinetic) to the gas specific energy as a function of radius for the 
70-80 lowest-mass clusters that have minimal substructure, in the 
112.9 /i~^Mpc simulations. The dashed, solid and dotted lines are 
for softenings of 10, 50 and 100/i^^kpc, respectively. 



Figure 3. The ratio of the gas specific energy (thermal plus ki- 
netic) to the dark matter specific energy within the virial radius 
for all clusters. In this and subsequent figures, clusters extracted 
from the 50, 112.9 and 153/i~^Mpc boxes are plotted using trian- 
gular, circular and square symbols, respectively. The 14 percent 
of the clusters with S > 0.2 are plotted with open symbols; the 
others with solid symbols. The best fit constant is show by the 
solid line. 



files show large variations in /3, we have averaged the profiles 
of the 70-80 smallest clusters, M < 2 x 10^" /i"^ M©, that 
show minimal substructure (5* < 0.2: see Section 2^) — these 
are the clusters with the lowest values of tr,min/io- There is 
no evidence for core heating in the runs with softenings of 
50 and 100/i~^kpc, whereas for the 10/i~^kpc softeing run 
/3 is suppressed not only in the core but throughout the 
cluster. This suggests that much of the heating is going on 
in sub-halos before the cluster forms. As the properties of 
these sub-halos are similar to those of the final cluster core, 
this would explain why the Thomas & Couchman formula 
(Equation ^) is more appropriate than that of Steinmetz & 
White (Equation]^. 

It is clear from Figures ^ and |2| that there is little differ- 
ence between the results for softenings of 50 and 100 /i~^kpc. 
They have similar values of /Jiso and the radial variation in 
f5 is minimal. We conclude that a softening of 50/i~^kpc is 
adequate for our purposes. 

The softenings for each of the other two boxes were 
chosen to give a similar central value of tr,min/io for the 
lowest-mass clusters as in the 50 /i~^kpc softening run. The 
values of /3i8o for each of the three production simulations 
are shown in Figure ^ Clusters with significant substruc- 
ture, S > 0.2, are shown as open symbols; smooth clus- 
ters with solid symbols. In this figure and all those that 
follow, clusters from the simulations of side 50.0, 112.9 and 
153.0 /i"^Mpc are shown as triangular, circular and square 
symbols, respectively. 

Navarro & White (1993) and Pearce, Thomas & Couch- 
man (1994) have shown that gas can gain energy from dark- 
matter during gravitational collapse even in the absence of 
numerical heating. The mechanism is quite straight-forward: 
both the gas and dark matter are stirred during the collapse, 
however the gas is able to thermalise its kinetic energy, then 
pick up more at the expense of the dark matter. There is 
no reason to suppose that this process should have occurred 
equally in low- and high-mass clusters. There are hints of 



a slight decline in (3iso at low masses in Figure |3|, but to a 
good approximation, /3i8o is independent of mass and equal 
to 0.94±0.03. Thus the fractional heating is comparable to, 
but slightly smaller than, the estimate of 10 per cent that 
we made from the Steinmetz & White formula. 

Because the gas has a greater specific energy than the 
dark matter, it is more extended and hence the baryon frac- 
tion within the virial radius is smaller than the global aver- 
age. We measure it to be approximately (0.85 ± 0.02) Qt in 
our simulated clusters, independent of mass. 



3 CLUSTER PROFILES 

In clusters, the dark matter is dynamically dominant. It 
is therefore important to have a good model of the dark- 
matter mass distribution. Navarro, Frenk & White (1995, 
1996, 1997) showed that the density profiles of simulated 
clusters in a wide variety of cosmological models are well- 
described by the formula 

^=^(iT^' 

where r = as is the radius and po and a are free parame- 
ters, and went so far as to describe this as a 'universal den- 
sity profile'. This formula must break down at large radii 
as it predicts infinite mass, but it appears to hold out to 
the virial radius, defined (in critical-density models) as the 
radius, riso, enclosing a mean overdensity of 180. 

In a previous paper Thomas et al. (1998) found that the 
NFW formula was indeed a good approximation to the mass 
distribution of an average cluster, but that there was a wide 
dispersion in the rate at which the density was declining at 
the virial radius. Hence we introduce a more general profile 



x{l + x)" 



(6) 
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Figure 4. The circular-speed profile for one of the clusters. The 
best-fit theoretical model is shown by the smooth curve. 



Figure 5. s versus Migo for all the clusters. The solid symbols 
represent clusters for which the slope of the density profile at two 
virial radii is within 25 per cent of the asymptotic slope at large 
radii. 



where s is a constant. Because the name is now so well- 
established in the literature, we call this the 'generalised 
NFW model, although we note that the case s = 3 was 
introduced first, by Hernquist (1990). 

The resolution of the simulations presented in this pa- 
per is not sufficient to determine the density profile in the 
cores of the clusters. This has been looked at in depth by 
Moore et al. (1998) and by Jing & Suto (2000). There is now 
evidence to suggest that the inner density cusp rises more 
steeply than p oc r~^, but that is of little consequence for 
the overall dynamics of the cluster and so we will stick with 
Equation B here. 



3.1 Best-fit circular speed profiles 

The density profiles of individual clusters are often very 
noisy which makes them difficult to match to any given theo- 
retical profile. It is better to use a cumulative profile such as 
the circular speed, Vc, which is anyway a dynamically more 
relevant quantity. For the above distribution. 



2 {l+xy-^ -{S-1)X~1 

Vr- oa — 



x{l + x) 



when s 7^ 2, and 



2 (1 -H3;)ln(l -hx) 



•(1 + ; 



(7) 



(8) 



when s = 2. 

An example of a circular-speed profile is given in Fig- 
ure ^together with the best-fit function of the form of Equa- 
tion [71 We have fitted and plotted the curve between twice 
the softening and twice the virial radius. This particular 
cluster was chosen because the goodness-of-fit as measured 
by the mean-square deviation from the theoretical curve is 
the median value for all the clusters. As can be seen, there 
are kinks in the profile, showing evidence of substructure: 
this is typical of most of the clusters in our catalogue. 

The best-fit values of s as a function of cluster mass 
are shown in Figure ^ Note that s has been limited to be 
less than or equal to 10. The asymptotic density profile at 



large radii has a slope of — (s -f 1). However, this is not 
always representative of the slope at two virial radii, the 
outer radius to which we fit the rotation curve. Hence we 
plot with solid symbols in the Figure only those clusters for 
which the two differ by less than 0.2s — this corresponds to 
a characteristic radius a < O.Srigo- 

A quick glance at Figure ^ makes the suggestion that 
s = 2 is a universal density profile seem surprising. However, 
it is hard to measure the density profiles in the outer parts 
of clusters with any degree of accuracy and the answer that 
one gets often depends upon the radial extent of the fit. At 
one virial radius, the slope of the density profile is far from 
— (s + 1) and so the asymptotic slope is poorly constained. 
In addition, s and a are strongly correlated and it is often 
possible to get a reasonable fit by forcing s = 2 and allowing 
a to vary. Hence the statement that the profile within one 
virial radius is 'consistent with and NFW profile' is largely 
meaningless. It is for this reason that we choose to fit the 
profile within two virial radii instead. 

We define clusters to be consistent with an NFW profile 
if the asymptotic slope of their density profile lies between 
—2.8 and —3.2 (i.e. 1.8 < s < 2.2). Just under a quarter of 
the clusters meet this criterion. 

The best-fit profiles of many clusters plotted with open 
symbols show a high value of s. However, this does not indi- 
cate steep density profiles at large radii because the best-fit 
core radii rise to compensate. Rather, it indicates that the 
functional form of the generalised NFW profile is a poor rep- 
resentation of the cluster. As an example consider the cluster 
shown in Figure ^. This is a smooth cluster: visually it ap- 
pears spherically-symmetric and its substructure parameter 
is small, S — 0.038. In addition, the best-fitting ellipsoid 
(see Thomas et al. 1998) is amongst the most spherical of 
any cluster with axial ratios of 1.14:1.0:0.94. However, the 
density profile, even out to 0.8 virial radii, is poorly fit by 
an NFW profile (dashed line). Also, it shows a sharp change 
in slope at this radius that cannot be matched by any gen- 
eralised NFW profile. The dotted line shows a generalised 
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Figure 6. The density profile for one of the clusters. The dashed 
line shows the best-fit NFW model. The dotted line shows the 
best-fit generalised NFW model with s = 10. The solid line shows 
a better, alternative model, as described in the text. 



Figure 7. The measured concentration assuming an NFW pro- 
file, xino, versus mass within the virial radius, MigOi for all the 
clusters. The 14 percent of the clusters with S > 0.2 are plotted 
with open symbols; the others with solid symbols. 



NFW model with s — 10. and a — S.lrigo- In order to re- 
produce the sharp decline in density at the virial radius, s 
has to be very large, but this then leads to a density profile 
that is declining far too rapidly in the outer parts of the 
cluster (and would steepen even more at radii larger than 
those shown in the Figure). A better representation of the 
density profile in this case is given by the dotted line which 
corresponds to the function 



Po 



:(1 +k2)s/2 



(9) 



where x = r/a as before, and s = 3.2, a = 0.40 rigo. This 
value of s is a much better estimate of the asymptotic slope 
of the density profile at large radii. 

It can be seen from Figure ^ that the line s = 3 roughly 
separates the solid from the open symbols in the upper half 
of the plot. Thus, where the outer slope of the density profile 
is well-defined, it generally lies between —3 (an NFW profile) 
and —4 (a Hernquist profile). The open symbols represent 
clusters, like that shown in Figure |^, that have a sharper 
break in their density profile than can be fit by a generalised 
NFW profile: these comprise about 30 per cent of the total 
cluster sample. 

The open symbols that correspond to values of s less 
than 2 are mostly clusters that show some degree of sub- 
structure, for which the density profile is not well-defined. 
These comprise another 16 percent of the cluster population. 

Although the spread in s is large, there is a weak trend 
for s to increase with mass. To make this more evident, 
we define an average low-mass and an average high-mass 
cluster by selecting all relatively smooth clusters, S < 0.2, 
in the mass ranges M < 1.2 x 10"/i~^ Moand M > 3.0 x 
10^'^h~^ Mq. The resulting profiles are extremely well-fit by 
our theoretical model with slopes of s = 2.0 and s — 2.3, 
respectively. 

It is clear from the above analysis that there is no 'uni- 
versal profile' for dark matter halos. A substantial propor- 
tion of clusters show obvious substructure, and even those 



that don't exhibit a wide variety of functional forms for the 
spherically-averaged density profiles of halos. 

3.2 Best-fit NFW profiles and concentrations 

The concept of a universal density profile is an attractive 
one. It makes modelling of observed clusters much sim- 
pler and it has the advantage that there is only one free 
parameter — the ratio of the characteristic radius in the 
NFW formula to the virial radius, xiso = rigo/a(s = 2).^ 
Therefore, we wish to see how well one can approximate 
cluster properties by assuming that they all foll ow the NFW 



3.1 



profile, in defiance of the results of Section 

The best-fit values of a;i8o as a function of mass are 
shown in Figure ^. There is a general trend of decreasing 
concentration as one moves to higher masses, but once again 
the scatter is large. The solid line in Figure ^ corresponds 
to the function 



2;i80 



3.7 



M 



1015/i-iM, 



(10) 



The 86 per cent of the clusters with S < 0.2 are spread 
equally above and below the line. We shall use this relation 
in the analysis that follows to see how well the simple NFW 
model predicts the measured scaling relations between tem- 
perature and mass. 

We can compare our concentration parameters to those 
of Navarro, Frenk & White (1997), by measuring halo mass 
in terms of M, , the mass of a spherical region with root- 
mean-square linear density fiuctuation of 1.69 today. For our 
simulations, M* « 3.3 x lO^^/i~^M0, and Migo/M, ranges 
from 3 to 300. This gives concentrations in good agreement 
with their CDM models. 

i NFW define X200 to be the 'concentration parameter', presum- 
ably using 200 as an approximation for the virial overdensity; in 
this paper we will use the term to stand for xigo instead — there 
is little difference between the two. 



© 2000 RAS, MNRAS 000, §-0 



A simulated rCDM cosmology cluster catalogue 7 




Figure 8. xigo versus S for all the clusters. Note that six points 
to the left and the bottom of this plot have been omitted for 
clarity. 
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Figure 10. Mean temperature profiles for (a) 79 low-mass, and 
(b) 34 high-mass clusters, as described in the text. The upper 
curve in each case shows the total specific energy of the gas and 
the lower curve the temperature, both in units of keV/k. 
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Figure 9. 5" versus Migo for all the clusters. A couple of points 
have been omitted from the bottom of the plot, for clarity. The 
solid line shows a median line that divides both the triangles and 
the circles plus squares into two equal groups. 



Interestingly, Figure |8| shows that the measured concen- 
trations are highly-correlated with the substructure statis- 
tic. Thus low concentrations do not normally occur in 
spherically-symmetric clusters, but come from spherical av- 
eraging of clusters that have substructure. This suggests that 
the weak trend of decreasing concentration with increasing 
mass is driven primarily by the fact that substructure in- 
creases with mass, as shown in Figure ^ 

We conclude this section by plotting in Figure |lo| the 
temperature profiles for the mean clusters described above. 
The lower curve in each panel shows the gas temperature in 
units of keV/fc, while the upper curve shows the total specific 
energy of the gas (thermal plus kinetic) in the same units. 
It should be noted that the gas is approximately isothermal 
within about one-third of the virial radius, but that its tem- 
perature drops rapidly in the outer parts of the cluster. This 



is a refiection of the fact that the specific energy of the dark 
matter declines in the outer regions of the cluster, because 
the ratio between the two is approximately constant out to 
two virial radii. 

It is a matter of some debate in the literature as to 
whether the observed temperature profiles of clusters are 
isothermal or decline to large radii. Using ASCA data, Ir- 
win & Bregman (2000) and White (2000) infer fiat or slowly 
increasing temperature profiles within 0.3ri8o, whereas 
Markevitch et al. (1998) find that the temperature drops 
by a factor of two within Q.Srigo- It is clear from our sim- 
ulations that isothermal profiles beyond about Q.3ri8o are 
inconsistent with our simple model and would probably re- 
quire some form of non-gravitational heating. 



4 PREDICTED SCALING RELATIONS 

In this section, we predict scaling relations for model clusters 
under the assumptions that they are spherically-symmetric, 
isolated and in hydrodynamical equilibrium. It is not that 
we are asserting that these assumptions are true, but rather 
we want to see how well such a simple model will perform 
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and to test the sensitivity of the scaling relations to changes 
in the model parameters. 

By scaling everything in terms of a characteristic ra- 
dius, a, and density po, for each cluster, we can make all 
the variables dimensionless. We write r = ax and p = pop, 
where 

and the tilde indicates a dimensionless quantity. Similarly, 
for the mass within radius r, 



m = 47rpo 1 
where 

rr /a 



m ■ 



px^dx. 



(12) 



(13) 



When we observe clusters, we usually measure their 
properties at or within the virial radius, r^ir = axA^ at 
which their enclosed density equals some multiple, Ac, of 
the critical density of the universe. For the rCDM model 
that we consider in this paper, Ac — 180, but we will set 
rs = axs and leave 5 as a free parameter, so as to be com- 
pletely general in our argument. 

From the above definitions, we can write po in terms of 

xs, 



2-3 

po = 60 (5i8o pc 

m6 



(14) 



where pc is the critical density, Siso 
rh{xs). 



5/180 and ms 



4.1 The virial temperature-mass relation 

We next calculate the velocity dispersion, ct as a function of 
radius for the NFW profile. To do this, we assume that the 
velocity ellipsoid is isotropic. It then becomes convenient to 
write = kT/nniH, where k is the Boltzmann constant, 
/im_ff is the mean particle mass for the gas, and T is the 
'dynamical temperature' of the halo. If the specific energy 
of the gas and dark-matter were the same, then T would be 
equal to the gas temperature. 

The spherically-symmetric Jeans' Equation for a fluid 
with an isotropic velocity dispersion is 



d(pr) _ mp 
dx x^ ' 

where 

k 



(15) 



(16) 



We are now in a position to calculate the mean dynam- 
ical temperature, Ts, within radius rs- By deflnition. 



rhsTs= / rdm= / pT x^dx. (17) 
Jo Jo 

Combining Equations |l^ and |l6| to eliminate a, and 
then substituting for po from Equation |l^, gives 



kTs 



GpruH (2407rpc <5i8o)^ ml 



16.0 5^30 



xsTs 
ms V1015/i-iMq 



keV, 



(18) 



E 




Figure 11. The shape factor xs Tx/ni/j in the dynamical 
temperature-mass relation of Equation ISl The dashed, solid and 
dotted lines correspond to s = 1.5, 2 and 3, respectively. 



where we have taken pmn ~ 10~^'^kg, corresponding to a 
fully-ionized cosmic mix of elements. 

When S — Ac, then Ts is the 'virial temperature'. We 
will use the notation Tigo rather than T^ir, to make explicit 
that we are referring to an overdensity of 180 (as in other 
cosmologies, the virial radius would correspond to a different 
overdensity) . 

If all clusters are self-similar (e.g. if they all follow the 
NFW model), and if we measure the total mass, ms, and 
mean dynamical temperature, Ts, within rs, then it is often 
assumed that they must obey the relation Ts cx m^^^ . This 
is not strictly true, however, because of the presence of the 
xsTs/rhs term which varies with concentration, xs- 

We show the variation of xs Ts/rhs with xs in Figure |ll| 
for values of s = 1.5, 2 and 3. We see that for values of a;i8o 
in the range 1.5 to 10 that are typical when fitting the NFW 
profile to clusters, a;i8o Tiso/jTiigo lies in the narrow range 
0.4-0.47. Fitting s = 1.5 profiles instead would give smaller 
values of a;i8o, and fitting s = 3 profiles larger values, so 
that these too give approximately the same answer. Thus 
the virial temperature-mass relation is robust to variations 
in the density profile. 



4.2 The gas density profile and the X-ray 
temperature-mass relation 

The virial temperature-mass relation predicted above is 
of little use observationally. Instead we need to use the 
emission-weighted. X-ray temperature. Previous studies 
such as Eke, Navarro & Frenk (1998) and Makino, Susaki 
& Suto (1998) have calculated the mean X-ray temperature 
and the luminosity of the intracluster medium assuming that 
it is isothermal and sitting in an NFW potential. Figure 
shows that the gas is far from isothermal and so we instead 
make the assumption, as indicated by Figure ^, that the 
ratio of the specific energies of the gas and dark matter is 
everywhere constant and equal to /3 = 0.94, i.e. Tgaa = T//3. 

We will further assume that the dark matter is dynam- 
ically dominant. Then 
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Figure 12. The shape factor x^Tx s/fhs in the X-ray 
temperature-mass relation. The dashed, soHd and dotted hnes 
correspond to s = 1.5, 2 and 3, respectively. 



1 



d(pgasT'gas) 



1 d(pT) 
P 



dx p dx ' 

which can be expanded to give 



d In pgas 
d In a; 



dlnp 
din a; 



(1-/3) 



dlnP 
d In X ' 



(19) 



(20) 



where P — pT. This shows that the gas is able to support 
itself more effectively and hence has a shallower density gra- 
dient than the dark matter. For /3 = 0.94, the asymptotic 
density gradient at large radii is —2.7, compared to —3 for 
the dark matter. 

The mean X-ray temperature within radius rs is defined 

by 



Lg Tx,s = 



where 



Tgas dL, 



L{x) = / pgas TgL dx 



(21) 



(22) 



and Ls = L{xs)- Note that we have assumed that the emis- 

~ 1 /2 

sivity of the gas scales as Tgas , although this is only really 
true at high temperatures where bremsstrahlung dominates. 
However, the dominant contribution to the variation in the 
integrand in Equation ^ comes from pgas and so this is a 
reasonable approximation. 

The values of the normalisation factor in the X-ray 
temperature-mass relation, xsTx,s/rhs, are shown in Fig- 
ure |l^. They can be seen to rise more steeply at large con- 
centrations than in Figure |l^. This is because the X-ray 
temperature declines far less steeply with radius than the 
dynamical one used previously. Figure |l^ suggests that the 
range of concentrations seen in our simulated clusters will 
add quite a lot of scatter to the mean X-ray temperatures 
within the virial radius, with the more concentrated clus- 
ters being hotter than the less-concentrated ones for a given 
mass by up to 40 per cent. 



Figure 13. Tigo versus Migo for a^U the clusters. In this and 
subsequent figures in this section, the 14 per cent of clusters that 
have a substructure statistic S > 0.2 are plotted with open sym- 
bols; others as filled symbols. The dashed line shows the best-fit 
power-law. The solid line shows the prediction from the isotropic 
NFW model. 



5 MEASURED SCALING RELATIONS 

5.1 The virial temperature- mass relation 

The location of each of our clusters on the Tiso-miso plane 
is shown in Figure Also shown by the dashed line is the 
best-fit power-law 



fcTiso ~ 6.9 



miso 



1015/1-1 M 



0.67±0.02 



keV, 



(23) 



which corresponds to a value of xsTs/rhg — 0.43. 

The predicted relation from Equation |l^, using the me- 
dian concentrations, xiso, taken from Equation |l^, is shown 
by the solid line in the Figure. It can be seen that this pre- 
dicts temperatures that are slightly too low, although the 
spread in temperatures is consistent with the spread in a;i8o 
seen in Figure ^. 

It would be possible to reconcile the predicted and mea- 
sured relations by using slightly larger values of a;i8o: Fig- 
ure H shows that the most regular clusters are biased to- 
wards higher concentrations. An alternative explanation for 
the low virial temperatures is that the dark matter velocity 
ellipsoid is not isotropic, as radial orbits are less efficient at 
supporting the particles within a given potential than trans- 
verse ones. This is discussed further in Section p. 3. 



5.2 X-ray versus virial temperatures 

The ratios of the dynamical temperatures to the X-ray tem- 
peratures, averaged within the virial radii, are shown in Fig- 
ure |l^. 

The prediction from the isotropic NFW model with con- 
centrations given by Equation ^ is shown by the solid line. 
It can be seen that the theoretical curve mimics the data in 
that it predicts that low-mass clusters should have lower val- 
ues of /3jf,i8o than high-mass ones. However the mean mea- 
sured values of /3x,i80 lie about 0.1 below the predictions. 
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Figure 14. The ratio of tlie dynamical temperature to the X-ray 
temperature within the virial radius for all clusters. The predic- 
tion from the isotropic NFW model is show by the solid line. 



Figure 15. The X-ray temperature versus mass within the virial 
radius. The dashed line shows the best-fit power law; the dotted 
line shows the best observational results from Horner, Mushotzky 
& Scharf (1999) and the solid line shows the prediction from the 
isotropic NFW model. 



Once again, this could be reconciled by using higher values 
of a;i8o or an anisotropic dark matter velocity ellipsoid. 

Although virial temperatures are notoriously difficult 
to measure, there is empirical evidence that the observed. 
X-ray temperatures of the intracluster medium are greater 
than the virial temperatures (e.g. Edge & Stewart 1991a,b; 
Bahcall & Lubin 1994). The latter give a value of /3x,i80 = 
0.94 for the highest-mass clusters, which is not too dissimilar 
from the values that we find in Figure However, the 
results of the next section suggest values that are lower than 
this. 



5.3 X-ray temperature-mass relations within a 
fixed overdensity 

The relationship between X-ray temperature and mass 
within the virial radius is shown in Figure The dashed 
line shows the best-fit relation from Evrard, Metzler & 
Navarro (1996, hereafter EMN), 



kTx 



7.85 



?ni8o 



1015/i-iMg 



keV. 



(24) 



They extracted 58 clusters in the temperature range l- 
10 keV, from three different sets of cosmological simulations 
with a variety of cosmological parameters. We have used the 
information given in their paper to interpolate their results 
to an overdensity of 180 — our results are in good agreement. 

The temperature normalisation is slightly higher than 
for the dynamical temperature-mass relation and corre- 
sponds to xsTs/rhs — 0.49. The prediction from Equation 
using values of xiso taken from Equation ^ is shown by the 
solid line on the Figure. Once again it gives temperatures 
that are slightly too low. 

The most extensive observational investigation of the 
Tx-Mi8o relation is by Horner, Mushotzky & Scharf (1999) 
who used many different ways to determine the mass within 
the virial radius. Their preferred measure, based on X-ray 
emissivity and temperature profiles, gives Tx oc M^gg , while 
those obtained from the isothermal /3-model and X-ray sur- 



face brightness deprojections are flatter, Tx oc Mj^gQ""'^®. 
They suggest that the slopes of these latter two measures are 
too low because they assume a dark matter density proflle 
of pdark oc whilst observations and simulations suggest 
that it is steeper, pdark oc r~^'*. 

The preferred relation from Horner, Mushotzky & 
Scharf (1999) is shown as a dotted line in Figure From 
this it is clear that either the measured cluster masses are 
too low, or more likely dissipationless simulations predict 
X-ray temperatures that are smaller than the observed val- 
ues. If we combine the results from Figures |l^ and ^ then 
this argues that observed high-mass clusters have values of 
Px,iso ~ 2/3. 

This high value of f3x, iso might be thought to pro- 
vide evidence for heating of the intracluster medium, but 
in fact radiative cooling can have the same effect! Pearce 
et al. (2000) have shown that the removal of low-entropy gas 
by cooling in the core of the cluster can raise the X-ray tem- 
perature of the cluster by 20 to 40 percent. This cooling is 
relatively more important in low-mass clusters and so would 
have the effect of flattening the Tx-M relation slightly. 

Mohr, Mathiesen & Evrard (1999) flnd in a sample of 45 
clusters that the ICM is more extended than the dark-matter 

1/2 

and that Tx 25 Afgas . In our simulations, we also flnd that 
the gas is more extended than the dark matter, due to its 
higher specific energy. However, the baryon fraction within 
the virial radius is approximately 0.85, independent of mass, 
and so obtain a temperature-gas mass relation that parallels 
the one for the total mass, Tx 2S Afgas^- Physical processes 
such as heating or radiative cooling are once again required 
(and act in the correct sense) to reconcile the observations 
and simulations. 

In Figure |l^, we show a similar plot to Figure |l^, but 
for the X-ray temperature-mass relation within an overden- 
sity contour of 1000. Once again, the isotropic NFW model 
slightly underpredicts the X-ray temperatures (but agrees 
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Figure 16. The X-ray temperature versus mass within a spheri- 
cal shell enclosing an overdensity of 1000. The dashed line shows 
the best-fit power law, the dotted line shows the observational 
results from Nevalainen, Markevitch & Forman (2000), and the 
solid line shows the prediction from the isotropic NFW model. 



Figure 17. X-ray temperature versus mass, for properties aver- 
aged within the Abell radius, for all the clusters. The dashed line 
shows the best-fit power law. The dotted line is the relation from 
Equation ^ for properties averaged within the virial radius. The 
solid line shows the corrected relation as described in the text. 



with the results of EMN). The best-fit power law is consis- 
tent with the self-similar prediction 

/ xO.67 

fcTx,iooo-13.0 (^^Qj^^^j keV. (25) 

The dotted line shows the observational results from 
Nevalainen, Markevitch & Forman (2000). The observed 
temperatures are again higher than the predictions from 
this non-radiative model but in a way that is now mass- 
dependent — this is consistent with our expectation that the 
effects of cooling and/or heating would be greater in lower- 
mass clusters. 



5.4 The X-ray temperature-virial radius relation 

From an observational point of view, it is much easier to 
measure the virial radius rather than the virial mass, and 
so it is surprising that more authors do not report the for- 
mer. There is no extra information here, however, as the 
two are simply related and so it is possible to rewrite Equa- 
tions ^ and ^ in terms of virial radius and X-ray temper- 
ature: 

(^) <»' 

™„«0.85(i|iji=!)*),-'Hi,c. (27) 

These follow with the self-similar form in that they have 
slopes of |. ri8o agrees with the results of EMN; riooo has 
a slightly lower normalisation. 

We note that the only observational measurement of 
''1000 by Vikhlinin, Forman & Jones (1999) is at first sight 
in gross disagreement with our prediction. They have 

rrooo « 0.415 (^^^) " h-'Upc (28) 



which gives a value at 5keV of 1.04/i~^Mpc — far in excess 
of the simulation results. The reason for this discrepancy 
is that Vikhlinin et al. have defined overdensity with re- 
spect to a mean baryon abundance of Q.bh^ = 0.010 which is 
much lower than either recent determination from primor- 
dial nucleosynthesis {Q,bh^ = 0.019, Tytler et al. 2000) or 
from X-ray measurements of the baryon fraction {Qth^ = 
0.0594/1^/^, Ettori & Fabian 1999). We can turn the argu- 
ment around and ask what baryon abundance would make 
our results compatible with Vikhlinin et al.. The radii they 
quote enclose a baryon overdensity of about 230 in our mod- 
els. Hence we require 

nth^ ^ 0.010 ^ 0.043, (29) 

This agrees with the value from Etorri & Fabian for h = 
0.54. 

5.5 Temperature-mass relations within the Abell 
radius 

If we measure cluster properties within a fixed radius, rather 
than a fixed overdensity, then the deviation from the self- 
similar scaling relations can be quite large. This is shown 
in Figure |l^ where we plot, for within the Abell radius, 
emission-weighted temperature versus mass. The best-fit 
power law, shown as the dashed line, has a slope of 0.81. 

The reason for the steeper slope is that, for low-mass 
clusters, the Abell radius is greater than the virial radius and 
so we are averaging properties over a larger volume than be- 
fore. This has the effect of lowering the X-ray temperature 
slightly (because the X-ray temperature is heavily weighted 
by emission from the centre of the cluster this effect is small) , 
but greatly increasing the mass. For high-mass clusters, how- 
ever, the virial radii are similar to the Abell radii and so 
there is no change. 

Because the X-ray temperatures are almost unaltered, 
the mass-temperature relation can be estimated by simply 
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Figure 18. X-ray temperature versus mass, for properties aver- 
aged within half the Abell radius, for all the clusters. The dashed 
line shows the best-fit power law. The dotted line is the relation 
from Equation ^ for properties averaged within the virial ra- 
dius. The solid line shows the corrected relation as described in 
the text. 



codes, yet yield very similar results. It should be noted that 
EMN do not quote results for an overdensity of 180, but 
they give sufficient information in their paper to allow one 
to interpolate from their preferred overdensity of 500. Once 
we do so, we find a Tx-M relation essentially identical to 
ours. 

Eke, Navarro & Frenk (1998) perform a similar study 
for 10 clusters in a A-CDM cosmology, ranging over a factor 
of only three in mass. Their results are consistent with the 
scaling relations but cover too small a dynamic range to 
provide a strong constraint. 

By contrast, Bryan & Norman (1998) use a smaller box 
and simulate a large number of clusters at relatively poor 
resolution. They find that /3 is an increasing function of clus- 
ter mass which suggests that numerical heating is having an 
effect. They also find a Tx-M relation that is flatter than 
the self-similar one, but only for an open cosmology. 

Yoshikawa, Jing & Suto (2000) have undertaken similar 
simulations to ourselves, but with a poorer spatial resolu- 
tion and covering a smaller mass-range. They too find the 
expected slope of I for the Tx-M relation. 



multiplying the right-hand-side of Equation |24| by the ratio 
of the mass within an Abell radius to that within the virial 
radius, MAbcii/A^iso. We have done this, using Equation 
to estimate the virial radii and assuming that the mass pro- 
file follows the NEW model. The result, shown as the solid 
line in Eigure |l^, closely follows the data points. 

Eigure |l^ is the same as for Figure ^ but for properties 
averaged within half an Abell radius. Once again the best-fit 
power law has a slope, 0.86, that is steeper than the self- 
similar relation, and the corrected NEW model provides an 
excellent fit to the data. 



6 DISCUSSION 

6.1 Comparison with previous simulations 

Previous authors have also looked at scaling relations in sim- 
ulations of clusters without radiative cooling. All agree that 
the gas is more extended than the dark matter, but disagree 
on the details of the Tx-M relation. The studies split into 
two types: those that simulate a small number of clusters 
in detail, and those that, like us, look at a large number of 
clusters at moderate resolution. 

Navarro, Erenk & White (1995) simulate only 6 clus- 
ters, extracted from a low-resolution simulation, but choose 
these to span a factor of thirty in mass and adjust the par- 
ticle mass so that each cluster contains several thousand 
particles within the virial radius. They find a Tx-M rela- 
tion in agreement with the virial relation and values of /3 of 
about unity (where their {3, as ours, includes a contribution 
from unthermalised motions in the gas) with no dependence 
upon mass. 

EMN have taken the Navarro, Frenk & White results 
and compared them with two other sets of simulations of 
isolated clusters from Metzler (1995) and Mohr et al. (1995). 
These span a range of cosmologies and use two distinct 



6.2 Beyond the isotropic NFW model 

It can be seen from Figures P and |l| that the isotropic 
NEW model underpredicts the virial temperature slightly 
and the X-ray temperature quite a lot. 

One possible explanation for the low virial temperatures 
is that the dark matter velocity ellipsoid is not isotropic, 
as radial orbits are less efficient at supporting the particles 
within a given potential than transverse ones. 

We have checked that a ratio of the transverse to radial 
velocity dispersion of at = 0.9 cjr is sufficient to bring the 
dashed and solid lines in Figures |l^ and |l^ into close agree- 
ment. We choose not to report on this in detail, however, be- 
cause we believe that it overly-complicates the simple model 
that we are trying to test. It should be remembered that few 
of our clusters are free from some form of substructure and 
this too would invalidate the simple NEW model. 

We prefer to conclude that the isotropic NEW model 
does pretty well, whilst bearing in mind that it underesti- 
mates the temperatures slightly. 



7 CONCLUSIONS 

We extract clusters, covering more than 2 decades in mass, 
from three simulations of the rCDM cosmology. This rep- 
resents the largest, uniform catalogue of simulated clusters 
ever produced. 

We fit the circular-speed profiles of the clusters with 
models corresponding to spherically-symmetric density pro- 
files of the form 



Po 



(30) 



^(l + a::)'' 

where s = 2 corresponds to the renowned NEW profile. 

We find that the scatter in the best-fit values of s is 
large, with under a quarter of the clusters being accurately 
fit by the NEW model. The others either have steeper outer 
slopes in their density profiles, show a sharp break in their 
density profiles that cannot be fit by the above form, or 
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have significant substructure. However, the mean, smooth, 
low-mass cluster does have s ~ 2, whereas the equivalent 
high-mass cluster has a steeper profile with s ~ 2.8. 

When we force s = 2, then the cluster concentrations 
show a large scatter, but the median concentration declines 
slightly with mass. This is driven primarily by an anti- 
correlation between concentration and substructure, with 
substructure being more prevalent in high-mass clusters. 

We investigate how well this median NFW model pre- 
dicts the cluster temperature-mass scaling relations, by 
deriving theoretical relationships between both dynamical 
and emission-weighted, X-ray temperature and mass in the 
isotropic NFW model, as a function of cluster concentra- 
tion. The virial temperature-mass relation, averaged within 
a spherical region enclosing an overdensity of 180, closely 
mimics the self-similar form and lies a few per cent above 
the NFW prediction. 

The gas in our clusters is hotter than the dark matter 
but we do not attach much significance to this: there will 
be a small amount of true and numerical heating due to 
heating by clumps and individual particles of dark matter. 
More importantly, we have chosen in this paper to neglect 
the effects of radiative cooling and heating associated with 
metal enrichment of the intracluster medium. 

When measured within spheres enclosing a fixed over- 
density, the X-ray temperature versus virial mass relation 
has a slope of approximately 0.67, in agreement with the 
self-similar prediction and with previous work (but extend- 
ing over a greater mass-range). The normalisation is lower 
than that of the observations. This is almost certainly be- 
cause we have chosen to neglect cooling and/or heating. 

The radius-temperature relation is in agreement with 
the measured value of Vikhlinin, Forman & Jones (1999) 
provided that the baryon fraction is large, as indicated by 
Ettori & Fabian (1999). 

When averaged within an Abell radius (or half an 
Abell radius), the X-ray temperature versus mass relation 
is steeper, with a slope of 0.81 (or 0.86). This is because 
the Abell radius is a greater multiple of the virial radius 
in low-mass clusters as compared to high-mass ones and so 
increases their measured mass. 

We have not commented on the X-ray luminosity in this 
paper, because our simulations are not able to fully resolve 
the X-ray emission in the cluster cores. This core gas has a 
cooling time that is anyway less than the age of the Universe 
and so cannot be correctly modelled using a non-radiative 
simulation. We note that the results of Pearce et al.(2000) 
suggest that radiative cooling will act so as to raise the 
temperature of the intracluster medium and so tend to bring 
our simulated clusters into agreement with the observations. 

In future papers we will contrast the dynamics of the 
gas and dark matter in clusters, consider the effect of heating 
and cooling processes, and compare the results from different 
cosmological models. 



ACKNOWLEDGMENTS 

The simulations described in this paper were carried out 
on the Cray-T3E at the Edinburgh Parallel Computing 
Centre as part of the Virgo Consortium investigations of 
cosmological structure formation. Interaction between au- 



thors was aided by a NATO Collaborative Research Grant, 
CRG 970081. OM is supported by a DPST Scholarship from 
the Thai government; PAT is a PPARC Lecturer Fellow; LO 
is a Daphne-Jackson Fellow, funded by the Royal Society. 



REFERENCES 

Bahcall N., Lubin L., 1994, ApJ, 426, 513 
Bryan G. L., Norman, M. L., 1998, ApJ, 495, 80 
Couchman H. M. P., Thomas P. A., Pearce F. R., 1995, MNRAS, 
452, 797 

Crone M. M., Evrard A. E., Richstone D. O., 1996, ApJ, 467, 489 
Edge A. C, Stewart G. C., 1991a, MNRAS 252, 414 
Edge A. C., Stewart G. C., 1991b, MNRAS 252, 428 
Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263 
Eke V. R., Cole S., Frenk C. S., Henry J. P., 1998, MNRAS, 298, 
1145 

Eke V. R., Navarro J. F., Frenk C. S., 1998, ApJ, 503, 569 
Ettori S. & Fabian A. C, 1999, MNRAS, 305, 834 
Evrard A. E., Metzler C. A., Navarro J. F., 1996, ApJ, 469, 494 
Governato F., Babul A., Quinn T., Tozzi P., Baugh C. M., Katz 

N., Lake G., 1999, MNRAS, 307, 949 
Gross M. A. K., Somerville R. S., Primack, J. R., Holtzman, J., 

Klypin, A., 1998, MNRAS, 301, 81 
Henry J. P., 1997, ApJ, 489, LI 
Hernquist L., 1990, ApJ, 356, 359 

Hjorth J., Oukbir J., van Kampen, E., 1998, MNRAS, 298, LI 
Horner D. J., Mushotzky R. F., Scharf C. A., 1999, ApJ, 520, 78 
Irwin J. A., Bregman J. N., 2000, ApJ, submitted (astro- 
ph/0003123) 

Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., Cole S., 
Evrard A. E., Yoshida N., 2000, MNRAS, submitted (astro- 
ph/0005620) 

Jing Y. P., Suto, Y., 2000, ApJ, 529, 69 

Lacey C. G., Cole S., 1993, MNRAS, 262, 627 

Lilje P. B., 1992, ApJ, 386, L33 

Makino N., Susaki S., Suto Y., 1998, ApJ, 497, 555 

Markevitch M., 1998, ApJ, 504, 27 

Markevitch M., Forman W. R., Sarazin C. L., Vikhlinin A., 1998, 
ApJ, 503, 17 

Metzler C. A., 1995, Ph.D. thesis, Univ. of Michigan 
Mohr J. J., Evrard A. E., Fabricant D. G., Gheller M. J., 1995, 
ApJ, 447, 8 

Mohr J. J., Mathiesen B., Evrard A. E., 1999, 517, 627 
Moore B., Governato F., Quinn T., Stadel J., Lake, G., 1998, 
ApJ, 499, 5 

Navarro J. F., Frenk C. S., White S. D. M., 1995, MNRAS, 275, 
720 

Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 563 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Navarro J. F., White S. D. M., 1993, MNRAS, 265, 271 
Nevalainen J., Markevitch M., Forman W., 2000, ApJ, 532, 694 
Pearce F. R., Couchman H. M. P., 1997, New Astron., 2, 411 
Pearce F. R., Thomas P. A., Couchman H. M. P., 1994, MNRAS, 
268, 953 

Pearce F. R., Thomas P. A., Couchman H. M. P., Edge A. C, 

2000, MNRAS, in press (astro-ph/9912013) 
Press W. H., Schechter P. G., 1974, ApJ, 187, 425 
Steinmetz M., White S. D. M., 1997, MNRAS, 288, 545 
Thomas P. A. et al.(the Virgo Consortium), 1998, MNRAS, 296, 

1061 

Thomas P. A., Couchman H. M. P., 1992, MNRAS, 257, 11 
Tytler D., O'Meara J. M., Suzuki N., Lubin D., 2000, PhST, 85, 
12 

Viana P. T. P., Liddle A. R., 1996, MNRAS, 281, 323 
Viana P. T. P., Liddle A. R., 1999, MNRAS, 303, 535 



© 2000 RAS, MNRAS 000, 



14 P. A. Thomas et al. 



Vikhlinin A, Forman W., Jones C, 1999, ApJ, 525, 47 
White D. A., 2000, MNRAS, submitted (astro-ph/9909467) 
Yoshikawa K., Jing Y. P., Suto Y., 2000 (astro-ph/0001076) 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^^X style file. 



© 2000 RAS, MNRAS 000, 



